Dynamic structure factor for '^He in two— dimensions 
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Recent neutron scattering experiments on "^He films have observed a zero-sound mode, its dis- 
persion relation and its merging with -and possibly emerging from- the particle-hole continuum.— 
Here we address the study of the excitations in the system via quantum Monte Carlo methods: 
we suggest a practical scheme to calculate imaginary time correlation functions for moderate-size 
fermionic systems. Combined with an efficient method for analytic continuation, this scheme affords 
an extremely convincing description of the experimental findings. 

PACS numbers: 67.30.ej, 67.30.em, 02.70.Ss 



Introduction: The two isotopes of Helium, '^He and 
''^He, give the opportunity to explore the quantum be- 
havior of many-body systems on a fundamental basis; at 
low temperature and pressures, they are the only neu- 
tral quantum liquids existing in Nature and an impres- 
sive complexity of physical phenomena is generated by 
merely pair interactions between particles and the effects 
of quantum statistics. In the investigation of the fasci- 
nating behavior of strongly correlated quantum systems 
a key role is naturally played by the low energy dynamics 
(see, for example, Ref. [l|). In addition, due to the very 
simple Hamiltonian, '^He and ^He many-body systems 
represent also extremely important reference models and 
test cases for general theoretical approaches.— 

Recently inelastic neutron scattering experiments have 
been performed on a monolayer of liquid "^He adsorbed 
on suitably preplated graphite: a collective zero-sound 
mode (ZSM) has been detected as a well defined excita- 
tion crossing and possibly reemerging from the particle- 
hole continuum typical of a Fermi fluidi^i^ From the the- 
oretical side, a quantitative description of such experi- 
mental findings has been achieved by a dynamical many 
body theory, without any adjustable parameters)^ The 
aim of this work is to undertake an ab-initio study of 
the ZSM in a strictly two-dimensional (2D) ■^He sample 
relying on Quantum Monte Carlo (QMC) methods. It 
has been shown that this ideal, strictly 2D model offers 
a realistic representation of the adsorbed liquid layer, as 
far as the liquid phase properties arc concerned.—"— 

The key quantity to be computed to compare with the 
ZSM observed in neutron scattering experiments on '^He 
systems is the coherent dynamic structure factor,— which, 
apart from kinematical factors, is related to the differen- 
tial cross section: 
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The brakets indicate a ground state or thermal average. 



H is the Hamiltonian operator, and — X^iLi e^*'' '"' is 
the local density in Fourier space. The ZSM of the system 
manifests itself in the shape of S{q,uj), appearing either 
as sharp peaks if it is long-lived or as broad structures if 
strong damping is present ii 



QMC methods give access to the coherent dynamic 
structure factor, S{q,uj), because they allow to eval- 
uate the intermediate scattering function F{q, t) = 

{e^^ pqe~'^^ P-^) by simulating the imaginary time dy- 
namics driven by the Hamiltonianj^iiS For a collection 
of ^He atoms, a very accurate microscopic description is 
afforded by the simple Hamiltonian 
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where is the mass of '^He atoms and v{r) is an effective 
pair potential among '^He atomsiii 

The correlation function F{q, r) is the Laplace trans- 
form of S{q,uj). Despite the well known difficulties re- 
lated to the inversion of the Laplace transform in ill- 
posed conditions, the evaluation of S{q,uj) starting from 
the QMC estimation of F{q, r) ^ has been proved to be 
fruitful for several bosonic systems J^ii^rJ^ 

For a Fermi liquid, the difficulty is further enhanced 
by the famous sign problem^ thereby the computational 
effort grows exponentially with the imaginary time and 
with the number of particles. While accurate approxima- 
tions exist to circumvent this problem in the calculation 
of static ground-state properties)^ we are aware of no ap- 
plications of approximate schemes such as the restricted 
patbi^ or constrained path^S methods to the calculation 
of imaginary-time correlation functions. 

Focusing on T = K, QMC calculations of ground 
state average replaces the unknown exact ground state 
V'o by the imaginary time projection of a trial function 

Vt,— '^'^ V'o = e~^"-ipT- The intermediate scattering 
function then reads: 
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Unfortunately, the projection time A required to filter out 
the exact ground state from the trial function is usually 
larger than the range of r needed to extract spectral in- 
formation, so that the total imaginary time 2A-|-r in Eq. 
([5]) is too large for practical purposes. In this paper we 



2 



- q = 0.267 A 



3 

S COS 




10 20 30 40 50 
0) (k) 




10 20 30 40 50 60 
CO (k) 



Figure 1: (Color online) Comparison between the spectral 
functions Si(g,a;) (filled curve down to the baseline) ob- 
tained with the GIFT algorithm from F\[q,T) and £'2(9,0;) 
(unfilled curve) obtained from F2{q,T) with GIFT for some 
wave-vectors q at density p — 0.047 A~^. The two spec- 
tral functions have a compatible shape, with a shift in energy 
compatible with — Eq . 



propose two related approximations which avoid the ex- 
tra time 2A, whereby the calculation becomes feasible for 
a few tens '^He atoms. The agreement with the measured 
dynamic structure factor is more than satisfactory. 

The dynamic fermionic correlation method: We adopt 
the standard Jastrow-Slater form for the trial function, 
= JD. The starting point of the present work is the 



following approximation: 
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where tpQ and ip^ are respectively the fermionic and 
bosonic ground states of the Hamiltonian, ([2]). Through- 
out this paper, a superscript F{B) indicates Fermi(Bosc) 
statistics and a subscript denotes the exact ground 
state. The convenience of the approximation ((4]) is that 
the extra projection time A does not compound the sign 
problem because it is applied only to the symmetric fac- 
tor of ipp- In the resulting approximate correlation func- 
tion 



Fi{q,T) = 
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the projection time between the determinants, which de- 
termines the severity of the sign problem, is limited to r. 
Fi{q,T) is an approximation of F{q,T), and its inverse 
Laplace transform, Si{q,uj), is an approximation of the 
coherent dynamic structure factor ([T]). The bias would 
vanish if "DipQ were the exact Fermi ground state. The 
QMC calculation of Fi requires the ratio of two correla- 
tion functions, Fi{q, r) = F2 {q, r) /FpQ (r), where 
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Both Fyq{t) and F2{q,T) are bosonic correlation func- 
tions and thus they can be evaluated with great accu- 
racy by means of exact bosonic QMC methods. -Fpq(t) 
is precisely the correlation function that was recently 
used in the fermionic correlations (FC) method^ to study 
the magnetic properties of ■^He films. On the other 
hand, F2 arises as an extension of the FC method 
to the calculation of the intermediate scattering func- 
tion, hence the name "dynamic fermionic correlation" 
(DFC) for the present methodology. Indeed, the func- 
tion F2 possesses very interesting features on its own: 
on one hand, it contains the exact fermionic spectrum, 
as can be seen from its spectral resolution: (<?,''■) = 

e^^(^" ^■^o where E![ are the fermionic en- 
ergy eigenvaules, Eq is the bosonic ground state en- 
ergy and 6„ = K/5_,-2? ^q^IV'^)] V(^o^|V'o^). If, more- 
over, the approximation dU is accurate enough, the co- 
efficients hn become, apart from an unessential normal- 
ization, the spectral weights of the exact intermediate 
scattering function, F{q,T). Therefore, computing the 
inverse Laplace transform of Fi and F2 we can obtain 
two different estimations for the coherent dynamic struc- 
ture factor, Si{q,uj) and S2{q,io + EQ —Eq) respectively, 
where the shift in cj comes from the definition of F2 in 
terms of the Bosc ground state. A robust test for the 
validity of the approximation ^ is at hand if it turns 
out that Si{q,uj) ~ S2{q,oj + Eq — Eq): as already 
noticed, F2 decays with the exact fermionic excitation 
energies (once shifted); moreover, if V^jjg has a small 
overlap on the fermionic excited states, it follows that 
g--rHj)^B quickly converge in r to e~'^^° . There- 
fore Fi{q,T) ~ e7^° F2{q,T), apart from an unessen- 
tial normalization. We have indeed verified that, in the 
present case, 5*1 (9, a;) and S2{q,oj + Eq — Eq) possess a 
very similar shape (see Fig[T]). 

Results: Wc studied a system of iV = 26 structure- 
less 1/2-spin fermions of mass m^, interacting via the 
Aziz potentialfii enclosed in a square box with periodic 
boundary conditions. Wc found in Ref. that this sys- 
tem size offers a good compromise between finite-size ef- 
fects and computational cost. Indeed the inverse Laplace 
transform becomes increasingly difficult as the range of 
fermionic energy eigenvalues relevant for the spectral re- 
construction departs from the reference energy of the un- 
derlying simulation, which is the bosonic ground state 
Eq-. this is precisely what happens as the system size 
increases because the gap Eq — Eq is an extensive quan- 
tity. The trial function tjj^ = J^D is the same as in 
Ref. 0, namely a two-body Jastrow factor J' and a Slater 
determinant D of plane waves with simple backflow cor- 
relations. We have focused on a density around 0.047 
A~^, close to the experimental conditions.— Moreover, 
we have explored the behavior of the sample at the den- 
sities 0.038 and 0.060 A~^ in order to investigate the 
density-dependence of the excitations of the system. In 
particular, the highest density was chosen very close to 
the freezing point The QMC evaluation of F2 requires 
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Figure 2: (Color online) From left to right the coherent dynamic structure factor, obtained as an average of several independently 
extracted Si{q,u!), for increasing wave vectors at p = 0.047 A~^. The orange shadows represents statistical uncertainties and 
filled (blue) circles are the available experimental data from Ref. and Ref. Q. The wave-vector shown in picture are those 
accessible from our simulation, the experimental wave vectors are q — 0.55 A~^(b), q = 1.15 A^^ (d), q — 1.25 A^^ (e) and 
q — 1.65 A~^ (f). We have used different scales in the panels to make more visible the comparison with experimental data. 
The (green) dashed line shows the dynamic structure factor of a fictitious system of bosons of mass ma. The bosonic peaks in 
the roton region are 5 — 9 times higher than the fermionic ones. 
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Figure 3: (Color online) Color map of normalized S2{q,i^) for many wave vector q. For better visibility, each 52 (g, a;) for 
different q has been normalized in order to have their maximum value equal to 1. The vertical scale have been shifted by a 
quantity — Eq , so that the excitation energies are measured with respect to the fermionic ground state. 



a simple generalization of the methodology we have fol- 
lowed in Ref. to compute Fp(j: a fictitious system 
of bosons of mass 7713 is simulated by an exact projec- 
tor Monte Carlo technique, the Shadow Path Integral 
Ground State.— The imaginary-time propagation was 
1.3125K^^ and the density matrix approximation was a 
Pair Product)^ with imaginary-time-step of 1/160 K~^. 
It is well known that, in order to extract information 
from imaginary time correlation function, an inversion of 
the Laplace transform in ill-posed contitions is necessary. 
This can be carried out by means of the Genetic Inver- 
sion via Falsification of Theories (GIFT) ^ which has 
been used to retrieve nontrivial spectral features in the 
study of low energy excitations of Bose superfluidai^'^'i^ 



and supcrsolidsji^ 

In Fig. [2] we show the comparison between our estima- 
tion of the dynamic structure factor of the '^He film and 
the experimental data.— The dynamic structure factor 
has been obtained as an average over several GIFT re- 
constructions of Si{q,uj) from independent estimates of 
Fi{q, r); this has made possible an estimation of the sta- 
tistical uncertainties which wc have shown in Fig. [2] with 
a yellow shadow. Wc note that the available experimental 
wave-vectors do not exactly match the reciprocal space 
grid defined by the simulation box. For q = 0.534 and 
q = 1.603 A^^, where the mismatch is minimal, a direct 
comparison is possible and the agreement is impressive. 
Inspection of the wave-vector dependence of the spec- 
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tra shows that the discrepancies seen at q = 1.069 and 
q = 1.336 are mostly due to the differences of q val- 
ues between theory and experiment. A major feature 
of the measured S{q,u;), captured also by the dynami- 
cal many body theory in Ref. 0, is the appearence of a 
low-energy peak for both small and large wavevector, in- 
terpreted in RefslH3 as a well defined collective mode, 
broadened in the intermediate q range because of mix- 
ing with the particle-hole continuum. In further agree- 
ment with the measurements, we find a similar behav- 
ior. Indeed the simulation can provide information even 
at small wave vectors, not accessible to the experimen- 
tal probe: at g = 0.267 the collective excitation (ZSM) is 
most pronounced, and the spectral weight of the particle- 
hole is negligible. It is remarkable that both the position 
and the shape of the calculated spectra have a physical 
meaning and are not artifacts of the reconstruction proce- 
dure. Further support to this conclusion is offered from 
a comparison with the dynamic structure factor of the 
fictitious ^He-mass bosonic system. The bosonic spec- 
trum has a completely different behavior, featuring an ex- 
tremely sharp peak with the usual phonon-maxon-roton 
dispersion, showing that the broadening of the fermionic 
spectrum is actually related to Fermi statistics. 

In Figl3] we report, in a color plot, the estimated 
52(9,1^). In agreement with the behavior of Si shown 
in Figl^l at low q we find well defined excitation energies; 
as the wave vector increases, again we observe broaden- 
ing of the ZSM. The ZSM dispersion E{q) as a function of 
the wave-vector q can be inferred from the distance be- 
tween the maximum in 5*2(9, w) and the value of the en- 
ergy gap between the Fermionic and the Bosonic ground 
state. The dispersion of the ZSM recalls the classical 
phonon-maxon-roton mode in bulk superfluid ^He; the 
roton energies decreases with the density and the min- 
imum moves to higher wave-vectors similarly to what 
happen for the bosonic liquid. Maxon energies increases 
with the density and also the zero sound velocity has the 
same trend. 

As a byproduct of the calculation of S {q,Lu) we also 
obtain the static density response function xiQ) from the 
moment -1, shown in Fig. We are unaware of previous 
QMC results for xil)- We have also calculated the static 
response of the fictitious bosonic system, which is signif- 
icantly less structured than the fermionic couterpart, as 
expected (see Fig. 

Conclusions: In this work we have presented an ab- 
initio estimation of the coherent dynamic structure fac- 
tor of 2D liquid ^He, a strongly interacting Fermi liquid. 



combining unbiased QMC sampling techniques with a 
statistical methodic for the analytical continuation from 
imaginary time to real frequencies under the only approx- 
imation dH) (Vj^ ~ T^iPn)- 

We find a well defined collective mode (the zero-sound 
mode) at small wavevector; its dispersion relation fol- 
lows a phonon-maxon-roton pattern, with a significant 
broadening in the intermediate wavevector range due the 
mixing with the particle-hole continuum. These fea- 
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Figure 4: The static response function of '^He obtained from 
Xq = -2p / dij (Circles) p = 0.038 A'^; (Squares) 

p = 0.047 k-^; (Triangles) p = 0.060 A'^. The filled circles 
show the results for a fictitious system of bosons of mass ms 
at p = 0.047 A-^ 

tures, including the shape and width of the spectra, 
are in close agreement with a recent neutron scattering 
experiment,—"^ and a sophisticated dynamical many body 
theory.— 

We expect that the approximation (jH) , as suggested by 
the accuracy of our results for 2D liquid '^He (See Fig.[T|), 
is a good one whenever the effects of quantum statistics 
are dominated by strong interactions. The formalism can 
be readily generalized by replacing the density fluctua- 
tions operator pg- to access other spectral properties such 
as spin fluctuations or particle-hole excitations. 
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